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1. Introduction 



Studies of tunneling dynamics in quantum models with triple-well potentials are becoming 
attractive in view of new trends in atomtronics associated with intensively discussed in 
literature [I], [2], j3], [I], (5], [6], [7], (8], (9) atomic diodes and transistors. Prospects to 
create new type devices are opened up with developing experimental tools for manipulations 



with cold atoms, and for laser monitoring and control 10 , 11 , 12 , 13 of tunneling 



transport. A triple- well structure has been used in experiments on light transfer in a 



triple-well optical waveguide 14 , 15 , classically modelling more subtle quantum coherent 
transport of neutral atoms or electrons in multi-well traps. 

Theoretical investments in atomtronic devices is often come to studies of tunneling trans- 
port in one-dimensional (ID) quantum-mechanical models, or models of Gross-Pitaevsky 
type, with double- or triple-well potentials. Properties of tunneling dynamics like particle 
transition between side wells with negligible filling of the central well, one way transport, a 



feedback of a small central well filling on the tunneling rate etc. 16 , 17 , 18 , 19 , 20 , 21 



are in the focus there. The resonance tunneling effect in the triple-well ID time-independent 



models and its experimental observations is discussed in [22], 23 as well. However, the 
analysis is carried out on exotic triple-well potentials, which are constructed by sewing 
separate (step-like or parabolic-type) wells or by separating wells with delta-function type 
barriers. As a consequence, spectra and wave functions in such potentials are recovered by 
numerical, i.e. non-analytical, methods. This makes difficult to study effects in the tunnel- 
ing dynamics associated with changes in the potential shape (for instance, with deformation 
of one well of a number of wells). 

Investigation of space-time evolution of quantum mechanical systems with potentials 
with several local minima is one of complicated and priority tasks. Complexity of the 
problem forces one to use various simplifications/approximations, e.g., matrix Hamiltoni- 
ans with account of few lowest levels [24], [25], [26], [27], or purely numerical simulations 
in phenomenological models. These simplifications can be used for figuring out common 
properties of tunneling dynamics, but not for studies of subtle effects. However, the de- 
tailed consideration of tunneling dynamics is becoming important in Condensed Matter 
Physics, where possibility of constructing new atomtronic devices based on cold atoms 
are intensively exploring. Therefore, studying subtleties of tunneling dynamics in multi- 
well potentials requires of using models with controllable characteristics, whilst studies of 
space-time dynamics in tunneling processes requires going beyond the standard approxi- 
mations (that take into account just few lowest states of the Hamiltonian and deal with 
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parameters chosen accordingly to the potential shape). The above mentioned requirements 
are inherent to exactly-solvable models of extended Supersymmetric Quantum Mechanics 



(SQM) with multi-well potentials 28 , 29 , 30 , 31 , 32 , 33 



Apparatus of extended N=4 SQM allows one to construct an exactly-solvable model with 
a multi-well potential, properties of which - eigenvalues of lowest states of the spectrum, 
degree of deformation of the potential etc. - can be preassigned. Technique developed 
in 34 , 35 can be used to derive the explicit analytical expression of the corresponding 



propagator. Hence, it is quite natural to extend methods of 32 , 33 to models with triple- 
well potentials. Free parameters of exactly-solvable models within N=4 extended SQM 
make possible to vary the shape of the potentials, hence the corresponding propagators, in 
a wide range, and in a consistent manner, studying effects of these changes on the dynamics 
of tunneling transitions. Recall that exact propagators take into account all states of N=4 
SQM Hamiltonians, that allows one to study tunneling dynamics of wave packages beyond 
the approximations of phenomenological approaches. 

The paper is organized as follows. In Section 2 we outline the procedure of construct- 
ing exactly-solvable isospectral Hamiltonians of N=4 SQM with three-well potentials, and 
discuss general properties of them. Previously obtained results of [33] are crucial to this 
end. In what follows we focus on harmonic oscillator model, and use the harmonic os- 
cillator Hamiltonian as a starting point in constructing a family of three-well potentials 
with different degrees of deformation. The corresponding exact propagators are obtained 
in Section 3. Gained results are applied in Section 4 to study features of tunneling dy- 
namics of the originally localized states in symmetric and asymmetric potentials. Systems 
with three-well potentials have a rich quantum dynamics, and we, in particular, observe a 
Josephson-type tunneling transition of a wave packet, the effect of its partial trapping, and 
a non-monotonic dependence of tunneling dynamics on the shape of a three-well poten- 
tial. Our conclusions are collected in Section 5, where we also discuss prospects of further 
development and applications of current results. 



2. Constructing isospectral Hamiltonians with triple-well potentials 

2.1. Warm up exercise: a double-well potential 

The procedure of constructing Hamiltonians with exactly-solvable N-well potentials from 
the initial Hamiltonian Ho, the spectrum and the wave functions of which are known, is 
based on adding complementary levels below the H ground state level with the energy E . 
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Mathematical basics of the procedure, which is known as the Crum-Krein's method, was 



proposed in original papers 36 , 37 , and intensively discussed in subsequent publications 



see, e.g., 34 , [35]). The method by Crum & Krein is equivalent to constructing isospectral 



Hamiltonians in polynomial SQM (reducible case) [28 31 , that is featured by successive 
application of supersymmetry transformations to the initial Hamiltonian. 

To construct Hamiltonians with triple-well potentials we will adopt previously obtained 



procedure 32 , 33 of getting exactly- solvable models of extended SQM with double-well 
potentials. Note, that the formalism of extended SQM is not principle to this end; the same 
results can be obtained by use of the Crum-Krein's method as well. However, some steps 
of the Crum-Krein procedure, such as, e.g., normalization of extra states wave functions, 
require a separate treatment with additional care. That is why we prefer to follow the 
method of successive application of supersymmetry transformations to the initial Hamilto- 
nian H Q , that directly leads to the isospectral Hamiltonians with multi-well potentials. 



Let us recall in brief main results of 32 , 33 , which will be in use in what follows. 



Suppose, we have an exactly-solvable (non-supersymmetric) model with Hamiltonian Hq = 
p 2 /2 + U(x) (in h = m = 1 units). For definiteness, let us fix H to be the harmonic 
oscillator Hamiltonian. In N=4 SQM, adding the extra level with the energy e below the 
ground state energy Eq of the initial Hamiltonian Hq, i.e. e < Eq, one may construct new 
isospectral Hamiltonian HZ , [j and its wave functions 



d 2 

H_ = (H -e)-—\ii(p(x,e,c), H -e = H + , 



I W ( (x, Ei) , tp(x, e, c) I 
i)- (x, Ei) = -. — ^ , , 2. 1 

N- 1 

i)_ (x, E = 0) = — r = 0{x, e, c). 

tp(x,e,c) 

Here ipZ.(x,Ei) are normalized wave functions that correspond to the states of the initial 
Hamiltonian H with energies E { (Ei ^ E ; energy is counted off e); W\^fi(x), f%(x) j = 
/1/2 — /1/2 i s the Wronskian of two functions /i^; /' = df/dx. ip and (p are new functions 
that appear in the following way. Solution to the auxiliary equation Ho(p(x,e) = e<p(x,e) 
for e < Eq consists of two linearly-independent non- negative functions ipj(x,e),I = 1,2 with 
the following asymptotic behavior: <pi(x — > —00) — > +oo,<£> 2 (x — > —00) — > 0, and (pi(x — > 
+00) — > 0,(p 2 (x — > +00) — > +00. Therefore, it is always possible to choose the integration 



1 Recall, that there are four Hamiltonians, Hf, H^, H + and H_ , in N=4 extended SQM, three of which, 
H+, and HZ, are independent (see, e.g., |32J for details). 
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constant c in such a way that the general solution ip(x,e,c) = N (<fi(x, e) + ccp 2 (x, £)) 
does not have any nodes on the whole axis. Taking the latter into account, the function 
tp(x,E,c) = N~ 1 /<p(x, e, c) is finite, and can be normalized to unity for every value of 
parameters e and c. Normalization of wave functions of excited states of HZ is preserved. 
For a specific choice of parameters e and c, the function (p(x,e,c), hence the potential of 
HZ, has local maxima and minima. The double- well structure of the potentials UZ{x,e,c) 
becomes sharp under (E — e)/E <C 1. 

also contains the Hi part, 



32 



33 



The super-Hamiltonian of N=4 SQM H£(a 1:2 = ±) 
that gets the extra level in the spectrum. Explicitly, the Hamiltonian H+ and its wave 
functions are 



Hi 




<p(x,e,c) 



l+A/*d££2(£,e,c 
1 



ipl{x,E = 0) 



N x ip(x,e,c) 



l + A Pd^ 2 (^,c)) 



d d 0(x,e,c) 
dx dx (1 + A d£ <^ 2 (£, £, c)) 



(2.2) 



Here A is the integrations constant, characterizing the general solution to the Riccati equa- 



tion (see Ref. 38 for details), whose value is restricted to A > —1. Note, that the po- 
tentials of HZ and H+ are shape invariant, i.e. UZ and are related to each other by 
non-coordinate transformations of their parameters. 



Fixing c = 1 in (|2.1|) and (|2.2|) and using the relation [39] 

( My\{x) - Aiy 2 (x) 



d£ 



W U(£),</ 2 (£) 



(A iyi (0 + A 2 y 2 (0Y 



{A{ + A\ 



\Aiyi(x) + A 2 y 2 (x) 



A 2 yi(xj) - Aiy 2 (xj) 
Amixi) + A 2 y 2 (xi) 
(2.3) 



we get, after inserting (2.3) in (2.1) and (2.2) 

d 2 



HZ = H~ 
A(e,X) = 

TpZ{x,Ei) 



dx 



2 In (tp 1 (x,e) + ip 2 (x,e)) 



h: 



dx 2 



A(oo,e,c = 1) - A - (A + l)A(-oo,e,c = 1) 
A(oo,e,c = 1) + A - (A + l)A(-oo,e,c= 1)' 

1 W[4jZ{x,Ei), v {x,e,l) 



In (tfi(x, e) + A(e, X)ip 2 (x, e)) , 
c<fx(x,e) - y 2 (x,e) 



A(x,e,c) 



y/2{Ei - e) <p(x,e,l) 

1 w(^Z{x,Ei),ip(x,e,A) 
sJ%Ei ~ e) ^,e,A) 



ipZ(x,E = 0) 



^X(x,E = 0) 



(pi(x,e) + cip 2 (x,e) 
N- 1 

= <p(x,e, 1) 



<p(x,e, 1) 

N A ' 
(p(x,e,A) 



<p(x,e,A). 
(2.4) 
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Clearly, the potential of HZ is defined by the symmetric combination ipi(x,e) + (p2(x,e), 
whilst U+ is defined by the asymmetric combination (pi(x,e) + A(e, \)<p2{x,e). The same 
concerns the wave functions of excited states of H~ and Hi. 



Eq. (2.3) is used to compute the normalization constants N 1 and A^ 1 of the corre- 
sponding ground state wave functions. They are related to each other by N^ 2 = AiV~ 2 = 
(A + l)N~ 2 , where 

N- 2 = --r- J 7T = W U(x,e ),ip2 (x,e )). (2.5) 

A(+oo, e, c = 1) — A(— oo, e, c = 1) V / 

For the initial Hamiltonian H of the harmonic oscillator model we get <£>i(£, e) = 

D v (y/2£), (p2(£,,e) = D v (—y2£), (£ = y/ujx,u = —l/2 + e,e = e/u; D v (x) is the parabolic 

cylinder function), and W^i,<^2j = 2y/mj/T(— v) with the gamma-function r(— v). The 

normalization constant N is 



A(+oo,f,i) r(-i/) 

Note that the only way of varying the shape of the potential in terms of dimensionless 
variables £ is to vary parameters A and e in the following range < e < 1/2, A>0. 
In natural units x, there is one more parameter u>, variation of which may change the 
potential shape, and, in particular, the positions of its local minima. 

2.2. Exactly -solvable Hamiltonians with triple-well potentials 

Now we are ready to proceed with constructing isospectral Hamiltonians with triple-well 



potentials. Eqs. (2.1), (2.2), (2.4) are basic equations to this end. 



Let's consider a double-well Hamiltonian as the initial one (cf.(2.4)) 



H -e 1 = H+= {H - e) - ^ In (x, e) + A(e, X)ip 2 {x, e)) . (2.6) 

The spectrum of H apparently contains the whole spectrum of the original Hamiltonian 
H , and the extra level with the energy e < E , which becomes the ground state of the 
Hamiltonian Hq. Clearly, Hq contains a deformed potential U+(x, e, A(e, A)). 

Next, consider solutions to the auxiliary equation H x( x ) = £ ix( x ) with S\ < e. As 
in the previous case of double-well Hamiltonians, one finds two linearly-independent (non- 
normalized) solutions to this equation 

W\<pi(x, e l ), <p(x, e, A) J WUp 2 (x, ei), <p(x, e, A) J 

X1(X ' £1 ' £ ' A) = j^A) ' ^,e 1)£ ,A) = , 

(2.7) 
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where ipi(x,£i),I = 1,2 are the solutions to equations H ipj(x,£i) = £i(pi(x,£i) with the 
original Hamiltonian H , and <p(x, e, A) = ^(x, e) + A(e, A)(y9 2 (x, e). 



Solutions (2.7) are non-negative, and have the following asymptotes: Xi( x ~~ >" — oo) — > 
+oo,X2(% — oo) and Xi( x ~~ ^ — > 0, x 2 (£ +oo) — > +oo. The minus 
sign in the expression of X2{x,£i,£,A) provides the correct asymptotic behavior of non- 
normalized functions. For the specific choice of integration constants, the general solution 

to H x(x,£i,£,A) = £ix{x,£i,£,A), £i < £, 

X(x, £%, e, A, c) = N(xi(x, £i, e, A) + cx 2 (x, £ U e, A)) 



does not have any nodes on the whole axis. 

In what follows we will set the constant c to unity. Then, the function 



x{x,e 1 ,e,A,c= 1) 



Xl(x,£i,£,A) + X2(x,£i,£,A) 



is finite, and it can be normalized to unity for every specific choice of the parameters. 
To construct a three- well Hamiltonian from the two- well Hamiltonian Hq, consider 



H Q -£ X = HX = H~(x,p) = X -p 2 + -[(V(ar,ei,e,A)) 2 +V'(x,£ 1 ,£,A) 



Here 



d d WU(x,ei,l),^ ( x,e,A 

V(x, £i, e, A) = — In x(x, £i, e, A, 1) = — In — ^ — 

dx dx (p{x,e,A) 



V= lx~( W ~l lnW ')' W(x) = -Un[l + \ I d£x 2 (£,£i,M) 



(2.8) 



with the superpotential W (see 38 for details), and (j)(x,£i,l) = <fii(x,£i) — <p 2 (x,£i)- 
The spectrum of HT completely coincides with that of Hf, and has the extra level with 



£ < Eq in compare to the spectrum of Hq. Recall, that the Hamiltonian H+ (2.4) has the 
potential with a double-well shape. 

The superpartner of HT is a triple-well Hamiltonian 



d 2 

H_ = H + - — In x(x, £i, £, A, 1) = (H - £ X ) 



dx" 1 



lnW ^ei.^^e.A) . (2.9) 



Its spectrum consists of the states of H , and two extra levels with the energies e and £\ 
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(Eq > e > Ei). The corresponding to these extra levels wave functions are 

N- 1 N- x ip{x,e,K) 



Xi(x,e 1 ,e,A)+X2(x,e 1 ,e,A) w(^>(x, e x , 1), <p(x, e, A)) ' 



(2.10a) 



2(e-ei) V rf2; x(^,^i^,A, 1) 



_ x^ 1 £i 1 £ 1 A 1 l) 
2 - ex) \dx x{x, £i, e, A, 1) 



W(0(x,£i,l),y?(x,£,A)^ 
(2.10b) 

ip+(x,Ei) 



Ei- e I s-sx ^[^+(x,Ei),ip(x,e,A) 
ip + (x,Ei) + — -0(x,5i, 1) 



Ei-Ex 



Ei — Ei 



W(0(x,ei,l),^(x,e,A) 



(2.10c) 



Here x{x,e u e,A, 1) = x(x, £i, e, A, c = 1) = iV(xi(a;, £i,£,A) + x 2 (ar,ei,e, A)). 



By use of (2.3) and the asymptotic behavior of Xi{ x i £ i5 e > A), / = 1,2 it is easy to 



calculate the normalization constant of the ground state wave function ^/q(x,Ei,e,A) 



N~ 



W(xi,X2j = 2(e-ei)W(^i(ar,ei),^ 2 (a:,ei)) . 



(2.11) 



To get a more symmetric representation of the extra levels wave functions, let's introduce 
another normalizations constant, N^ 2 — > 2(e — ei)N^ 2 . Then eq. (2.10b) becomes 

AT A V(x,£i,l) 



#l(x,£,£l, A) 



W(0(x,£i,l),^(x,5,A) 



Remind, solutions tpi(x,ei),I = 1,2 enter (|2~8), (2.9), (2.10) through £i, 1) = 
</?i(x,£i) — (p2(x,ei) = 0(x,£i,Ai = 1). Below we will show that fixing the starting Hamil- 
tonian H to the harmonic oscillator one, a symmetric potential UZ of HZ corresponds to 
the choice A = Ai = 1. Deformations of UZ , even if A x = 1, mainly affect side wells, and 
to lesser extent the central well. 



2. 3. Shape-invariance of triple-well potentials in N=4 SQM isospectral Hamiltonians 
Let us establish the shape-invariance of the potentials entering HZ and H+ (see 



32 



33 



for a double-well case). H_ and H + = H — £\ are related via (cf. (2.2)) 



" f//n l} ^'H' • A./;;>/a^(c;,,A.. : 



(2.12) 
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Using (2.3), one gets 



x(x,£ X ,£,A,l) 



1 + A ff d£ x 2 (£, ei, e, A, 1) Xi(x, £ X ,£, A) + (A + l)x 2 (s, £ X , e, A) 



, N A 2 = (A + 1)A> 



As a result 



d 2 



d 2 



HX = H-- — \yi \ X (x, e x , e, A 1? A) | = (H -e x )- — In W(0(x, e ls A x ), <^(x, e, A) J , (2.13) 
where Ai = A + 1 (A > —1 for normalized wave functions), and x{ x -> £ i> e i Ai, A) = 



iV(xi(x, £i, e, A) + A 1 x 2 (a ; , £i, £, A). Comparing eq. (2.13) with eq. (2.9), one concludes 
on the shape-invariance of HZ and Hf. 



Wave functions of H+ states comes from the H_ states wave functions (2.10), replacing 
(x, ei, 1) with <f>(x, Ei, Ai) = (p x (x, Ex) - Aiip 2 (x, e x ) 



*o(^,£i,Ai;£,A) 
V x (x,e,k;e x ,kx) 



Nll<p{x,E,h) 



WU(x,£ X ,A x ),(f(x,£,A) 

iVAV(x,ei,Ai) 
W(0(x,ei,A 1 ),^(a;,£,A) 



(2.14a) 
(2.14b) 



£/j — E\ 



Ei-e I e-e x W[il> + (x,Ei),<p(x,e,A) 
ip + (x,Ei) + — —<f>(x,e x ,Ai 



Et — Ei 



W[(f){x,E x ,A X ),(p(x,£,A) 



(2.14c) 



Two remarks are in order now. Eqs. (2.9), (2.10) and (2.13), (2.14) include the wave 



functions of H , as well as non-normalized solutions to Schrodinger equation H ip(x) = 
£ip(x) with e < E . And potentials entering HZ and H+ depend on a number of parameters 
(e x , £, A x , A, . . . ) that makes it possible to vary their shapes, from purely symmetric to 
essentially deformed. 

2.4- Triple-well model from the harmonic oscillator Hamiltonian 
Consider the initial Hamiltonian to be the harmonic oscillator one 



p 2 u 2 x 2 



2 2 

As it has been noted above, non-normalized solutions of Hn<f(x) = £ip(x),£ < En are the 
parabolic cylinder functions 



the Wronskian of which is 

Then, the normalization constants are 



N~ 



4(1/ 



'71LU 



and 



AT 



4(z/ - /i) 



r(-/i) 



A^r 2 = AiV -2 , A>0, z/ = -^+e 



A^ 2 = AiN~ 2 , Al > 0, = + £! 



£1 



Fig. (JlJ) shows the shape of the potential Uz(£, fA, v. A) for change in the energy of the 
ground state for a fully symmetric case (see the left plot on Fig. ([!])), and for a non-trivial 
value of the deformation parameter of the side wells (see the right plot on Fig. ([I])). For 
a symmetric potential, increasing the value of S\ enlarges the width of the central local 
minimum, and reduces the widths of the side local minima. The distance between the 
side wells gets also increased. Hence, the energy e\ is the parameter which controls the 
shape of potential. Note, that varying e± does not affect the energies of the first and of 
the second excited states. Having a non-trivial deformation (A < 1), decreasing e\ results 
in the growth of the barrier in front of the more deeper minimum. 



U 



u 





4 

a) b) 
Fig. 1: UZ(X,H, u,A) shape variation for: a) A = 1, v = -0.02; b) A = 0.05, u = -0.02. 



It is important to trace back the behavior of wave functions corresponding to lowest 
states of HZ (under-barrier states). When all wells have the same depth (see the left panel 
on Fig. ([2])), one may notice that the ground state wave function (solid curve) is almost 
entirely concentrated in the central well, and is completely absent in the side wells. Whilst 
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the wave functions of the first excited state (dashed curve) and of the second excited state 
(dotted curve) are almost zero in the central well. Their values in the under-barrier region 
is enough to provide the tunneling process. In the case of A ^ 1 the wave functions 
distribution does not essentially change in the central well. A redistribution of the first 
and of the second excited states wave functions turns out in the side wells: the right well 
essentially contains the first excited level wave function, while the second excited level wave 
function is mostly concentrated in the left well (see the right panel on Fig.p])). 




Fig. 2: U_ (£,/!, is, A) and wave functions for: a) A = 1, /i = —1, v = —0.02; b) A = 0.05, 
H=-l,v= -0.02. 

A more detailed analysis of HZ and H+ wave functions behavior is postponed to Section 

4. 



3. Quantum-mechanical propagator in systems with three-well potentials 

Study of temporal evolution in quantum-mechanical systems with multi-well potentials is 
one of the priority tasks of Quantum Mechanics. Due to complexity of tunneling dynamics 
in real systems one has to use different simplifications or approximations, such as two- 
mode approximation, or phenomenological matrix models based on account of a few levels. 
These models can be used in studying general properties of tunneling dynamics, leaving 
apart more subtle quantum-mechanical effects. 

A most general way to describe space-time evolution in quantum-mechanical systems 
is to use the quantum-mechanical propagator. Recall that the exact propagator contains 
contributions of all states of the quantum-mechanical Hamiltonian. Suppose, a quantum- 
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mechanical object, say, a wave packet with the squeeze parameter R, 



$(x ,0) 



we 



2R\ 1/4 



7T 



exp 



1 



/ \2 

[X — x ) e 



(3.1) 



is initially, i.e. at t = 0, localized at xq. Then, its space-time evolution is described by |42| 

/+oo oo 
rfx A(x, t; so, 0)$(x , 0), A(x, t; x , 0) = £ ^ n (a?X(ar )e- <i3D ' . (3.2) 
00 n=0 

Here K(s,t;x o ,0) is the propagator, the explicit form of which allows one to study the 
dynamics of localized states in potentials of any complexity. Currently, there are only 
few examples of exactly-solvable models with explicitly derived expressions of propagators 
K(x,t;xo, 0) [43]. As a rule, they correspond to Hamiltonians with one-well potentials. 
A procedure of getting new models with propagators related to that of exactly-solvable 



models with one-well potentials was proposed in 34 , 35 . Following 34 , 35 , the exact 
expressions of propagators for quantum-mechanical systems with two-wells potentials were 
derived by two of us in [33], and the results were applied for studying the tunneling 
dynamics of localized states. Here we extend the method of 1 33 1 to Hamiltonians with 
triple-well potentials. 

Let's introduce the propagators KZ and corresponding to Hamiltonians HZ (2.9) 
and H+ (2.13). The explicit form of can be obtained on account of the shape-invariance 
established in the above. Having expression of K+, KZ is derived by setting Ai = 1 in the 
expression of K+. The latter is 

K+(x,t;x o ,0) = J2K( x > E n)K( x o,E n )e-* Eot 

n 

+*i(s, e, A; £l , Ai)*i(ar , e, A; e u A^e - * 6 * + V (x, e x , An e, A)^ (^o, e u A l5 e, A)e~ i£lt , 

(3.3) 

where the summation (in the absence of the continuous spectrum) is carried over the states, 
related to the initial Hamiltonian Hq. Other terms of (3.3) correspond to two extra levels 



with energies e and E\ [e > ex). To establish the relation of the propagator (3.3) to that 
of the exactly-solvable model with the Hamiltonian Hq we will only consider the first term 
on the l.h.s. of (3.3). Wave functions of H+ and H Q are related to each other with 

1 



2 y /(E n -e){E n -e 1 
where L x and L x are the N=4 SQM intertwining operators 

Lx = lx~ llx^^' 6 ^^ 



L x L x i/j + (x,E n ), E n y£e,ei 



(3.4) 



11 



d d 



d d 



+ — In V (x, £1, A x ; e, A) = — + — In 



cp(x,e,A) 



dx dx 



dx dx 



W( ei, Ai), cp(x, e, A) 



By use of (3.4) the first term on the l.h.s. of (3.3) is written down as 

1 



J2 K( x > E n )^X(x , E n )e- iEnt = L X L X0 L X L X0 

n n 

On account of the following identity 

1 1 



±{E n -e){E n -e x ) 



ip + (x, E n )ip + (x , E n )e 
(3.5) 



-iE n t 



(E n — e)(E n — ei) (e — £1) \E n — e E n — e\ 
and a property of a propagator K (x, t; z, 0) 



/+00 
dzK o (x,t;z,0)t/;-(z,E n ) 
-00 



one gets 



V *+(rc, £ n )*+(x , £ n )e~^ = x Xo x . Xo / dz K (x, t; z, 0) (G (z, x ; e) - G (z, x ; e,)) 

„ 4 ( £ - £ l) ./-do 



(3.6) 



where 



G Cz,a;o;e) = 22 



ip + (z, E n )if) + (x , E n ) 



E n — e 

is the Green function of H . It is well known (see, e.g., [44]), the Green function admits 
the following representation in terms of solutions to Schrodinger equation on interval [a, b] 

2 



Gjx,y,e) 



W /,,/, 



e)/ r (3/, e)0(3/ - x) + e)/ r (x, e)0(x - y)) , (3.7) 



satisfying the boundary conditions fi(a,e) = and f r (b,s) = 0. In our conventions 
fi(x,e) = (f2(x,e), f r (x,s) = <fi(x,s)] the same is for e\. 



Substituting (3.7) into (3.6) and acting with L Xo L XQ , after some tedious calculations we 
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get 

K( x > E n)K( x o> E n )e~ iEnt = L X L X (FxO, x , e, A, e u A x ) + F 2 (x, x , e, A, e ls A x )) 

n 

Fxfoa^e.A.eijAx) = - 1 rK (x, t; x , 0) ? 0(^0, £iAi) x 



x / dzK (x,t; z, 0) (A0 2 (z, e)6>(a;o - z) - (p 1 (z,£)8(z - x )) 

J — oo 

F 2 (x, x , e, A, ei, Ai) = — — rAT (x, So, 0) 7 y?(x ,£iAi) 



x / dz K (x, t; z, 0) (Aiy? 2 (z, ei)0(x o - z) + y>i(z, ei)0(z - ar )) • 

(3.8) 



Note, that every term on the l.h.s. of (3.8) contains Ko(x,t;xo,0), but with the opposite 
signs. As a result, these terms compensate each other, and we avoid the need to operate 
(differentiate) with a singular function (distribution). The multiplier in front of the integral 
in Fi (F 2 ) coincides, modulo the normalization constant, with the wave function of the first 
excited (of the ground state) level of H+. Therefore, the final expression of the propagator 
has the following form 

^(x, t ; Io ,0) = -M ! ^fkMi iiiX 

r+00 

x / dz K (x,t;z,0) (A<f 2 (z,e)9(xo — z) - ipx(z,e)9(z - x )) 



#o(zo,ei, Ai;e, A) = 



K 1 



x / dzK o (x,t;z,0) {A 1 (p 2 {z > e 1 )0{x o - z) + <p 1 (z,£i)9(z - x )) 
J —00 

+ ^x(x, e, A; e u A 1 )^ 1 (x , e, A; e u A 1 )e~ iEt + V (x, £1, A l5 e, A)V (x , £1, A x ; e, A)e~ ieit . 

(3.9) 

Choosing the initial exactly-solvable Hamiltonian Hq we have to specify the solutions 
<fj(x,e),I = 1,2, and the propagator K (x, t; x , 0). For the harmonic oscillator model 
tpi(£,e) = D u (y/2£), <p 2 (£,,e) = D v (-\/2£), and K o (x,t;x o ,0) is as follows 



43 



we -Mr( s +n)\ / 4 - w \ 

— — : exp — Ux + x ) cosut - 2xx \ , (3.10) 

27rsma;T / \2smut L J J 
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where t = nix/u + r, n G N, 0<r< n/u. The presence of exp(— Z7r(| + n)) in (3.10) 
provides the correct behavior of the propagator for all values of t. 



4. Dynamics of localized states in multi-well potentials 

Now we are ready to study features of the localized states dynamics in three-well potentials, 
with taking into account the contribution of all states forming a localized state $(x, 0). 



Basic equations are eqs. (2.14), (3.2), (3.9). Changing the parameters (e, £i,A, Ai) enables 
us to investigate the effect of the potential shape on the tunneling dynamics in a wide range. 
To specify the consideration we will fix the initial Hamiltonian to the Hamiltonian of the 
harmonic oscillator model, in terms of dimensionless variable £ = \fujx. Clearly, changing 
the value of u results in the appropriate change of the potential shape (in particular, in the 
localization of the side minima) in natural units x. Furthermore, it is interesting to study 
features of the tunneling dynamics in dependence on the choice of $(a;,0). This becomes 
especially true for amplification of tunneling transitions by use of a small part of the initial 
packet in the central well. 

In general, dynamics of localized states cannot be correctly described with taking into 
account just few lowest under-barrier states. The contribution of higher excited states 
in forming of $(x,0) is becoming essential during the tunneling process, and their role 
increases with the degree of localization. The approach we follow enables to consider the 
features of the localized states dynamics with taking into account all states of an exactly- 
solvable Hamiltonian with symmetric and asymmetric multi-well potentials. 



To proceed further, let us write down the basic equations (2.14), (3.2), (3.9) in terms 
of dimensionless units £. We have 

/+oo 
-oo 
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where 0) is the wave packet at the initial time, and 



J — oo 



#o(fo,ei,Ai;e,A) f 



#o(£, r ; c, 0) (Ai/y-v^O^o - C) + D^V2Q6{C - U 

-oo 

+ , e, A; e l9 AO^x^o, e, A; e 1 ,A 1 )e~ i£T + £i, A i; e, A)^ (£o, £i, A x ; e, A)e" 



■ieiT 

(4.2) 



with 



ue -iTr{\+n) \ 1/2 



exp 



[(e 2 + C 2 )cosT-2eC] ; (4.3) 



27rsinr / r V2sinT 
T = nn + r, n <EN, < t < n. Wave functions of the ground and of the first excited states 



(2.14a), (2.14b) are, in terms of dimensionless units, 

A^>(£,£,A) 



¥ (f,ei,Ai;e,A) 

*i(^e,A;ei,Ai) 

where = A^" 2 , N^ 2 = AA>- 2 , 
N~ 2 = 4(z/-/i) : 



W(0(£,£i, AO, A) 
Wto, £1^0,^,6, A) 



(4.4a) 
(4.4b) 



7TW 



AT" 2 =4(l/-^): 



7TW 



A, Ai > 0, n — —1/2 + Si/u, v = —1/2 + e/u. Functions (ft and entering (4.4), are as 
follows 

ei, A x ) = D M (V20 " A 1 D /1 (-V20, e. A ) = D v (y/2g) + AD U (-V2£) . 



Now eqs. (4.1), (4.2), (4.3) and (4.4) become basic equations to study the features of the 



tunneling dynamics in a wide range of varied parameters (e, E\, A, Ai, u). 

We get started with the initially localized in one of side wells wave packet 

,2H\ V4 



$(C,0)=^J exp(-i(C-Co)V^) 
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where Co is the position of the corresponding well minimum. Next, we are going to figure 
out the effect of the deformation degree on the tunneling dynamics, and its dependence 
on energies of extra levels in A, £ 1; Ai). Looking at Fig.([lJ one may notice that 

the three- well structure is most evident for the choice of parameters e, £j, when the state 
with the energy e forms the tunnel doublet with the initial Hamiltonian Hq ground state 
energy E , A = (E — e)/uj <C 1, and the energy S\ is fixed to be |ei| ~ uj. Therefore, 
it is advisable to fix the value of A, and to study effects of changing E\ on the tunneling 
dynamics. Decreasing the ground state energy of results in deepening the central 
well, and in concentration of the ground state wave function there. However, the ground 
state wave function of still plays an essential role in the tunneling dynamics, since 
transition rates from side wells depend on the overlap integrals between the ground state 
wave function and wave functions of other states. In this respect, changes in E\ with 
fixed values of A enable to change the tunneling rate. What concerns deformations of the 
potential, we will mainly consider them by use of A variations. It leads mostly to changes 
in the shape of side wells, though a deformation of the central well is also takes place (see 
Fig. p| 



4-1- Symmetric case 

All states of Hamiltonians HZ and H+ form the complete and orthonormal basis. We 
would like to show the efficiency of the basis in series expansion of $(£, 0) over the lowest 
states in symmetric and asymmetric potentials (states counting starts from the ground 



state), comparing the exact results of (4.1) with the following approximation 

Umax 

mz,t)\= cn^: ) ) +(e,K)ex P (-^ n T) 

n=0 



(4.5) 



The squeeze parameter value R (cf. (3.1)) is fixed by the requirement of minimal number 
of states giving the essential contribution to the series expansion of $(£,0). 

Table ([TJ contains the coefficients of the $(£, 0) series expansion over the states of 
H(£,n, v, A = 1) for different choices of (fi,u,R). "Left" and "Right" indicate the minima 
in which the wave packet is initially localized. 

It is important to note that once the energies of extra states are close to the ground 
state energy of the harmonic oscillator model, the main contribution to the series expan- 
sion of $(£, 0) comes from three lowest states. This result is in favor of application of 
phenomenological three-level Hamiltonians to tunneling processes at fixed parameters of 
the potential 24 , 25 , 26 , 27 . Considered in Table (fl| values of parameters correspond 
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state number 


1 


2 


3 


4 


5 


6 


7 


8 


9 


A — 1,R— 1.0, Left 
H = -0.03, v = -0.02 


0.56 


0.699 


-0.417 


0.038 


-0.048 


0.06 


-0.067 


0.062 


-0.042 


A — 1,R — 1.0, Left 
H = -0.032, v = -0.02 


0.54 


0.698 


-0.441 


0.041 


-0.052 


0.066 


-0.072 


0.066 


-0.045 


A—1,R— 0.66, Left 
/i = = -0.02 


0.038 


0.68 


-0.69 


0.162 


-0.138 


0.098 


-0.051 


0.016 


-0.003 


A = l,i? = 0.65, Right 
= -2.013, i/ = -0.02 


0.014 


-0.69 


-0.69 


-0.175 


-0.145 


-0.094 


-0.037 


0.007 


0.008 



Table 1 



to a wide central well (see the left plot on Fig.Q), and to the strong overlapping of wave 
functions, that leads to the appearance of a sufficient portion of the wave packet in the 
central well in the tunneling process (see Fig.® 



1*1 



A 



\ 



0S>s 




a) 



b) 



Fig. 3: |$(£, T)| dependence in the symmetric case: Plot a) corresponds to the 1st row of 
Table (HI); Plot b) corresponds to the 2nd row of Table (fT|. 



It is worth mentioning that a small change in the \i value (the second row of Table ([T]| 
results in essential rearrangement of dynamics of tunnel transitions. This is indicated by 
increasing the revival time of the wave packet, as well as by more complicated structure of 
|$(£,T)|. Therefore, we conclude on non-monotonic dependence of the tunneling dynamics 
under changing properties of a three-well potential. It means, that the values of parameters 
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of phenomenological models 24 , 25 , 26 , 27 should also have a non-monotonic behavior. 



When the ground state energy |ei| ~ u, the contribution of the ground state wave 
function in the series expansion of $(£, 0) is very small (cf. 3rd and 4th rows of Table (pj)). 
Indeed, looking at Fig. ([2]), one may notice that the ground state wave function is located 
at the central well, and its value in the side wells is negligibly small. Space-time evolution 
of the wave packet is characterized by a very small portion of |$(0,T)| in the case, and by 
independence of the revival time on the choice of the ground state energy. In other words, 
a transport of "particles" from side to side wells, with negligibly small central well filling, is 
realized in the case. The transport time is 1/2 T rev i va i = ir/A, and the complete returning 
time of the packet is equal to T reviva i, due to time reversibility of quantum mechanics. To 
avoid the wave packet returning, laser beams action on a quantum-mechanical system was 
proposed in 16 -18 . Clearly, this mechanism breaks the time invariance. Note also, that 



the central well filling rate decrease under increasing Hence, changes of enable to 
control the rate of the central well filling in the tunneling processes. 




a) 




b) 



Fig. 4: |$(£, T)| dependence in the symmetric case: Plot a) corresponds to the 3rd row of 
Table (HI); Plot b) corresponds to the 4th row of Table (HI). 



4-2. Asymmetric case 

Let's consider the effect of deformation of the potential on the tunneling dynamics of wave 
packets. We will suppose Ai = 1, and will vary the potential shape by changing A. As it 
has been noted above (see Fig.Q), the most significant changes in the shape is achieved at 
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A — > 0. Table ([2]) contains the coefficients of the series expansion of 0) over the states 
of (£, /i, v, Ai = 1,A = 0.05). Note, that a good approximation of $(£, 0) is achieved 



CTQTO Till T VI r~\ DT' 


i 
i 


o 


Q 
O 


A 


o 


u 


7 


o 


Q 


A = 0.05, # = 0.6, Left 
// — _o 03 77 — —0 02 

Lai U.UUj f U.U/j 


0.768 


0.21 


-0.564 


0.034 


-0.026 


0.016 


-0.006 


0.005 


-0.022 


A = 0.05, i2 = 0.8, Right 
H = -0.03, v = -0.02 


0.17 


-0.96 


-0.128 


-0.015 


-0.014 


-0.011 


-0.007 


0.009 


-0.021 


A = 0.05, R = 0.3, Left 
(i = -l,u= -0.02 


0.078 


0.209 


-0.938 


0.136 


-0.072 


0.053 


-0.074 


0.11 


-0.12 


A = 0.05, R = 0.5, Right 
[i = = -0.02 


0.012 


-0.95 


-0.215 


-0.044 


-0.029 


-0.012 


-0.005 


-0.02 


0.059 



Table 2 



by taking into account ten lowest states of H+. Other states contribute no more than 5% 
of probability density of the initial wave packet for different values of parameters, entering 
the potential and the wave packet. Hence, we conclude that the proposed basis is relevant 
to describe the dynamics of localized states even in the case of a strong deformation of a 
three-well potential. 



Space-time evolution, described by eqs. (4.1), (4.2), (4.3), (4.4) for a deformed potential, 
gets a number of differences in compare to the case of a symmetric potential. First of all, 
the tunneling transitions dynamics is characterized by partial trapping of the wave packet 
at the initial well (see Fig.(|5])). This effect consists in suppressing tunneling transitions in 
another side well, and it is strongly pronounced when the wave packet is initially situated 
in the more deep well. The mechanism of the partial "trapping" of the wave packet in the 
initial well is quite simple. If the maximal contribution to $(£,0) located at one of the 
side wells comes from the under-barrier state, the value of its wave function in the other 
side well is very small. Hence, the contribution of this state to the tunneling to the other 
side well process is insignificant. Other under-barrier states give small input into tunneling 
transitions, due to their minor contributions in forming $(£, 0). 

What concerns the central well filling, in the case of the initial location of $(£, 0) in 
the left well (see the left panel on Fig.Q) its value is big, while in the case of location 
of $(£, 0) in the right well, the central well filling value is small (see the right panel on 
Fig. (O). The same situation is observed for \n\ ~ 1 (see Fig. (pi)), though the portion of the 



wave packet at the central well is negligibly small even when the wave packet is initially 
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a) 



b) 



Fig. 5: |$(£, T)| dependence in the asymmetric case: Plot a) corresponds to the 1st row 
of Table p|; Plot b) corresponds to the 2nd row of Table pi). 



located at the left well (compare with the previously considered left panel on Fig.rt5 



1*1 



200 



300 




a) 




b) 



Fig. 6: |$(£,T)| dependence in the asymmetric case: Plot a) corresponds to the 3rd row 
of Table pi); Plot b) corresponds to the 4th row of Table pi). 



As for the case of a symmetric potential, the dynamics of tunneling transitions strongly 
depends on small changes of the tunnel doublet value Ai = 1/2 — e/u = —v. In particular, 
it leads to increasing the period of the wave packet revival time T revival = 2tiuj/{Eq — e). 
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5. Conclusions 



In this paper we have developed the previously proposed approach [33], based on exactly- 
solvable quantum-mechanical models with multi-well potentials and the corresponding exact 
propagators, to study the dynamics of initially localized states in triple-well potentials. 
By use of the procedure of constructing isospectral N=4 SQM Hamiltonians with multi- 
well potentials, it has been obtained general expressions of Hamiltonians with three-well 
potentials, and their wave functions. Parameters entering the obtained expressions enable to 
change the shapes of three-well potentials in a wide region. We have got general expressions 
of quantum-mechanical propagators, which take into account all states of the Hamiltonians. 

Specifying the initial Hamiltonian to that of the harmonic oscillator model, we applied 
the results to study features of the tunneling dynamics. It is important to note the relevance 
of the basis which has been used in the series expansion of the initial wave packet 0) 
to this end. This basis is formed by states of the obtained Hamiltonians under arbitrary 
deformation of the potentials. Even in the case of a large squeeze degree of the wave 
packet, taking ten lowest states is sufficient to approximate T) in a correct way. This 
claim has been verified by comparing (4.5) to computed with the exact propagator (4.1). 

Tunneling dynamics in a symmetric three-well potential demonstrates a number of in- 
teresting properties. When the tunnel doublet states, which are formed by the first and the 
second excited levels, are far from the ground state, the wave packet T) demonstrates 
a "Josephson"-type behavior. The revival time period T reviva i is determined by the value 
of the tunnel doublet A = (Eq — e)/oj in the case. Furthermore, the central well filling is 
small during the tunneling transport from side to side wells, and the filling increases with 
growth of |ei|. Therefore, changing the ground state level one may control the central well 
filling during the tunneling transport. 

When three lowest states are close, the central well filling becomes large. The wave 
packet revival time T reviva i, as well as its qualitative behavior, strongly depends on changes 
of Ai = e — S\. Non-regularity in the behavior of T reviva i is explained by its dependence 
as on A" 1 , as well as on Aj" 1 . The value of T reviva [ is defined by common multiple of A -1 
and A^ 1 . Hence, small deviations in one of them (say, in Ai) results in the significant 
increasing T reviva i, that leads to the essential changes in |<3>(£, T)| behavior. 

Tunneling dynamics of wave packets in the deformed three- well potentials demonstrates, 



as in the case of two-well potentials 33 , the effect of partial trapping of the wave packet 
in the original well. Put it differently, the portion of the wave packet tunneling from one 
to another side wells is sufficiently small, independently on the energies of lowest states. 
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The central well filling is essential only when the lowest states are close to each other, and 
when the wave packet is initially located at a less deep well. Tunneling in the deformed 
potentials possesses a non-regularity of |$(£, T)| too. 

Studies of tunneling in atomtronic devices are essentially based on a time-dependent 
interaction, provided, for instance, by the laser-beam radiation. Effects of having such 



interactions were investigated in systems with two-well potentials 45 , 46 , 47 , 48 , 49 



where it was demonstrated the efficiency of the tunneling dynamics control by use of 
external driving forces. For systems with three-well potentials theoretical modeling the 
effects of driving forces are at an early stage, and are based on matrix Hamiltonian models 



27 , 49 , or on potentials, properties of which can be studied by computer simulations 20 



Handling exactly-solvable models with three-well potentials with widely varied shapes, and 
with the corresponding exact propagators, open up the possibility to study the tunneling 
dynamics in time-dependent external forces. We believe that the obtained results will be 
useful in filling a gap in these investigations. 
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